Dependence of boundary lubrication on the misfit angle between the sliding surfaces 
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Using molecular dynamics based on Langevin equations with a coordinate- and velocity-dependent 
damping coefficient, we study the frictional properties of a thin layer of "soft" lubricant (where the 
interaction within the lubricant is weaker than the lubricant-substrate interaction) confined between 
two solids. At low driving velocities the system demonstrates stick-slip motion. The lubricant may 
or may not be melted during sliding, thus exhibiting either the "liquid sliding" (LS) or the "layer 
over layer sliding" (LoLS) regimes. The LoLS regime mainly operates at low sliding velocities. We 
investigate the dependence of friction properties on the misfit angle between the sliding surfaces and 
calculate the distribution of static frictional thresholds for a contact of poly crystalline surfaces. 

PACS numbers; 46.55.-|-d, 81.40.Pq, 61.72.Hh, 68.35.Af 



I. INTRODUCTION 

The problem of boundary lubrication is very interest- 
ing from the physical point of view and important for 
practical applications, but it is not fully understood yet 
[l|, 01 . Conventional lubricants belong to the type of liq- 
uid ("soft") lubricants, where the amplitude of molecular 
interactions within the lubricant, Vu, is smaller than the 
lubricant-substrate interaction, Vsi. Due to strong cou- 
pling with the substrates, lubricant monolayers cover the 
surfaces, and protect them from wear. A thin lubricant 
film, when its thickness is lower than about six molecu- 
lar layers, typically solidifies even if the conditions (tem- 
perature and pressure) are those corresponding to the 
bulk liquid state. As a result, the static friction force 
is nonzero, fs > 0, and the system exhibits stick-slip 
motion, when the top substrate is driven through an at- 
tached spring (which also may model the slider elastic- 
ity). In detail, at the beginning of motion the spring 
elongates, the driving force increases till it reaches the 
static threshold fs- Then a fast sliding event takes place, 
the spring relaxes, the surfaces stick again, and the whole 
process repeats itself. This stick-slip regime occurs at low 
driving velocities, while at high velocities it turns into 
smooth sliding. 

Since the pioneering work by Thompson and Rob- 
bins [3, 4] , who studied the lubricated system by molec- 
ular dynamics (MD), the stick-slip is associated with 
the melting-freezing mechanism: the lubricant film melts 
during slip and solidifies again at stick. Such a sliding 
may be named the "liquid sliding" (LS) regime. However, 
at low velocities the "layer over layer sliding" (LoLS) 
regime sometimes occurs, where the lubricant keeps well 
ordered layered structure, and the sliding occurs between 
these layers 0. 
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In real systems the substrates are often made of the 
same material and may even slide along the same crys- 
tallographic face, but can hardly be assumed to be per- 
fectly aligned, especially if the substrates have polycrys- 
talline structure. In the majority of MD simulations, 
however, both substrates are modelled identically, i.e., 
they have the same structure and are perfectly aligned. 
This fact may affect strongly the simulation results, as 
became clear after predicting the so-called superlubric- 
ity, or structural lubricity [5]. For example, the "dry" 
contact (no lubricant) of two incommensurate rigid infi- 
nite surfaces, produces null static friction, = 5-^. 
If the surfaces are deformable, an analog of the Aubry 
transition should occur with the change of stiffness of 
the substrates (or the change of load Q): the surfaces 
are locked together for a weak stiffness, and slide freely 
over each other for sufficiently high stiffness (this effect 
was observed in simulation [6|). 

In a real- world 3D contact, incommensurability can 
occur even for two identical surfaces, if the 2D sur- 
faces are rotated with respect to each other. Simula- 
tions [lol - [T^ do show a large variation of friction with 
relative orientation of the two bare substrates. Sim- 
ilarly to the ID Frenkel-Kontorova system, where the 
amplitude of the Peierls-Nabarro barrier is a nonana- 
lytic function of the misfit parameter, in the 2D system 
the static frictional force should be a nonanalytic func- 
tion of the misfit angle between the two substrates. This 
was pointed out by Gyalog and Thomas in their study 
of the 2D FK-Tomlinson model [l5|. However, surface 
irregularities as well as fluctuations of atomic positions 
at nonzero temperature makes this dependence smooth 
and less pronounced. For example, MD simulations fl^ 
of the Ni(100)/Ni(100) interface at T = 300 K showed 
that for the case of perfectly smooth surfaces, a 7r/4 ro- 
tation leads to a decrease in static friction by a factor of 
34-^330. However, if one of the surfaces is roughened with 
an amplitude 0.8 A, this factor reduces to 4 only, which is 
close to values observed experimentally. Miiser and Rob- 
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bins [6| noted that for a contact of atomically smooth 
and chemically passivated surfaces, realistic values of the 
stiffness usually exceed the Aubry threshold, thus one 
should expect fs = for such a contact. An approxi- 
mately null static frictional force was indeed observed ex- 
perimentally in the contact of tungsten and silicon crys- 
tals More recently the friction- force microscopy ex- 
periment made by Dienwiebel et al. [l^ demonstrated 
a strong dependence of the friction force on the rota- 
tion angle for a tungsten tip with an attached graphene 
flake sliding over a graphite surface, where sliding occurs 
between the graphene layers as relative rotation makes 
them incommensurate. 

The case of lubricated friction was investigated by He 
and Robbins [l^, [l^l for a very thin lubricant film (one 
monolayer or less). The dependences of the static [l^ 
and kinetic (20j friction on the rotation angle were calcu- 
lated. The authors considered the rigid substrates of fee 
crystal with the (111) surface and rotated the top sub- 
strate from (/) = to 7r/6. It was found that static friction 
exhibits a peak at the commensurate angle (0 = 0) and 
then is approximately constant; the peak/plateau ratio 
is about 7 (for the monolayer lubricant film, where the 
variation is the strongest). The kinetic friction varies 
slowly with a minimum at the commensurate angle and 
a smooth maximum aX 4> ~ 7r/18 = 10°, changing by a 
factor near two. Also, the kinetic friction decreases with 
velocity at (/> = 0, while it increases at the other angles. 

The goal of our work is a detailed MD study of stick- 
slip and smooth sliding for lubricated system with ro- 
tated surfaces. Compared to the work by He and Rob- 
bins 1^ 2^, we study thicker lubricant films, up to five 
atomic layers thick. We explore a fairly basic model (see 
Sec, im and Ref. (2l| . Interactions are of simple Lennard- 
Jones type, typically each substrate consists of two layers 
with 12x11 atoms, arranged as a square lattice, and the 
lubricant has 80 atoms per layer) which allows us a rather 
detailed study of the system dynamics for long simula- 
tion times. This model attempts to address the effects 
of relative crystal rotations in generic lubricated sliding, 
without focusing on a specific system. While the mi- 
croscopic interactions are treated at a minimal level of 
sophistication, we describe energy dissipation by means 
of a "realistic" damping scheme, with a damping coeffi- 
cient in Langevin equation, which mimics the energy ex- 
change between the lubricant atoms and the substrate. 
This is certainly important for smooth sliding, the ki- 
netics of melting and freezing processes and the ensuing 
metastable configurations which emerge at stick during 
stick-slip regime. 

Section IIIII presents typical simulation results. The 
model exhibits stick-slip at a low driving speed, which 
changes to smooth sliding with increasing speed. In the 
smooth sliding regime, as well as during slips in the stick- 
slip regime, the system indeed exhibits either the LS or 
LoLS regime, depending on the simulation parameters, 
and in particular the rotation angle. A new important 
result of the present study is that the LoLS regime should 



be observed much more often than the LS regime. Sec- 
tion |IV] discusses and summarizes the obtained results. 



II. THE MODEL 

As we address rather general properties of lubricated 
friction, we explore a relatively simple model in simula- 
tion. This allows us to span a wider range of sliding ve- 
locities and longer simulation times as well as to analyze 
the atomic trajectories in greater detail than by simu- 
lating fully realistic force fields appropriate to a specific 
system. The MD model was introduced and described 
in previous work [3. [2l|: therefore, we only discuss here 
its main features briefly. We study a film composed by 
few atomic layers and confined between two substrates, 
"bottom" and "top" . As illustrated in Fig. [TJ each sub- 
strate consists of two layers. The external one is rigid, 
while the dynamics of the atoms belonging to the inner 
layer, the one directly in contact with the lubricant, is 
fully included in the model. The rigid part of the bot- 
tom substrate is held fixed, while the top substrate is 
mobile in the three space directions x, ?/, z. 

All atoms interact with pairwise Lennard- Jones poten- 
tials 



V{r) = Vaa' [raa'/r] - 2 (raa'/r) 



(1) 



where a,a' = s or I for the substrate or lubricant atoms 
respectively, so that the interaction parameters Vaa' and 
faa' depend on the type of atoms. Between two sub- 
strate atoms we use Vss = 3 and an equilibrium distance 
rss — 3. The interaction between the substrate and the 
lubricant is much weaker, Vsi = 1/3. For the "soft" lu- 
bricant itself, we consider Vu = 1/9 and an equilibrium 
spacing ru = 4.14, which is poorly commensurate with 
r^s . The equilibrium distance between the substrate and 
the lubricant is rgi = ^{rss + Hi)- For the long-range 
tail of all potentials we adopt a standard truncation to 
r < r* = 1.49 r/;. The atomic masses are m; — rus — 1. 
The two substrates are pressed together by a loading 
force // per substrate atom (typically we used the value 
// = 0.1). All parameters are given in dimensionless units 
defined in Ref. [2l| , for example the model units for force 
is Vss/tss- The chosen parameters correspond roughly to 
a typical system where energy is measured in electron- 
volts and distances in angstroms, so that forces are in 
the nanonewton range. 

The main difference between our technique and other 
simulations of confined systems lies in the dissipative cou- 
pling with the heat bath, representing the bulk of the sub- 
strates. We use Langevin dynamics with a position- and 
velocity-dependent damping coefficient ?y(di, Vi), which is 
designed to mimic realistic dissipation, as discussed in 
Refs. [3, HJ]. In a driven system the energy pumped 
into the system must be removed from it. In reality 
energy losses occurs through the excitation of degrees 
of freedom not included in the calculation, namely en- 
ergy transfer into the bulk of the substrates. To model 
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Figure 1: (Color online) Three typical stick configurations of the A^; = 3 system during stick-slip motion for misfit angles 
(a) (j) = 0°, (b) 17.7°, and (c) 28.7°. Lubricant atoms are dark/red, deformable substrate atoms are clearer/blue, and fixed 
substrate atoms are white. Driving speed is = 0.03, load /; = 0.1. 



this fact, damping should occur primarily when a moving 
atom comes at a small distance di from either substrate. 
Moreover, the efficiency of the energy transfer depends 
on the velocity Vi of the atom because Vi affects the fre- 
quencies of the motions that it excites within the sub- 
strates. The damping is written as r](d,v) = rii(d) r]2{v) 
with 771 (d) = 1 — tanh[((i — d*)/d*], where d* is a charac- 
teristic distance of the order of one lattice spacing. The 
expression of 772 (w) is deduced from the results known 
for the damping of vibrations of an atom adsorbed on a 
crystal surface (see Ref. and references therein). 

In simulations we explore the "spring" algorithm, 
where a spring is attached to the rigid top layer, and 
the spring end is driven at a constant velocity. We ap- 
ply periodic boundary conditions (PBC) in the x and 
y directions. The geometric construction of the rotated 
substrate is explained in the Appendix □ In simulation, it 
is simpler to rotate the bottom substrate only. The initial 
configuration of the lubricant is prepared as a set of Ni 
{Ni — 1, 2, 3, and 5) closely packed atomic layers. Most 
simulations include 80 atoms in each lubricant layer, al- 
though we increased the system size by up to 16 times 
to check for size effects. The system is then annealed, 
i.e., the temperature is raised adiabatically to T ~ 0.6, 
which exceeds the melting temperature Tm {T„i ~ 0.1 
for our lubricant ^2^) and then decreased back to the 
desired value. After preparation of the annealed configu- 
ration, we perform a standard protocol of runs: starting 
from the high-speed v = 1 LS regime, corresponding to 
a sliding speed comparable to the sound velocity of the 
lubricant in its solid state, we reduce the driving velocity 

V in steps down to the value v = 0.01 (which produces 
stick-slip motion for most employed simulation param- 
eters), and then increase w, up to v = 3. Typical MD 
snapshots are shown in Fig. [TJ 

To estimate fs in the stick-slip regime we select the 

V = 0.03 runs and take an average of the peak spring 
force immediately prior to slip. A lower sliding speed 
would lead to slightly larger friction, due to longer ag- 
ing of pinning contacts, but would require longer simula- 
tion times to record the same number of stick-slip events. 
The moderate speed v = 0.03 realizes a fair compromise, 
which in practical simulation times produces a ^ 10% 
underestimate of fs with respect to its value at adiabat- 



ically slow sliding. To find the kinetic friction force 
in the smooth sliding regime, we average the spring force 
over the whole run. 



III. SIMULATION RESULTS 

The simulation results are summarized in Fig. [2j Qual- 
itatively, the results agree with those of He and Robbins 
(19I [20I , with stick-slip motion at low driving velocities 
and smooth sliding for u > 0.1. At zero substrate tem- 
perature, the static friction can vary with the misfit angle 
<f> by two orders of magnitude, and the kinetic friction for 
smooth sliding at low velocity (e.g. as shown for v = 0.3 
in Fig. [2Jd) by more than one order of magnitude. For the 
thin lubricant films, Ni < 3, the static friction peaks for 
perfectly aligned substrates, ^ = 0. This does not occur 
any more for thicker films. The friction achieves sharp 
minima at the angles (j) = 17.7°, 25.7° and 40.5° as will 
be discussed below. At large driving velocities, v > 1 
(which is in fact huge, comparable to the solid-lubricant 
sound velocity) the lubricant film is completely molten, 
and friction becomes almost independent of 0. 

The variation of friction with is the most pronounced 
for the one-layer lubricant film. The simulation results 
for this thinnest film are presented in Fig. |3l The static 
and low-speed {v < 0.3) kinetic friction force displays 
sharp minima for the angles (f> = 17.7°, 25.7° and 40.5°. 
For these "special" angles, the lubricant film remains or- 
dered and slides together either with the top or the bot- 
tom substrate both during slips and in smooth sliding (we 
call this regime as the "solid sliding", or SS). Of course 
the motion is not rigid but corresponds to a "solitonic" 
sliding mechanism |24l - [26| . For the other angles stud- 
ied, at stick configurations the film orders locally, with a 
structure adjusted partly to the bottom and partly to the 
top substrates, while during slips, as well as at smooth 
sliding, the film is 2D- melted (LS regime). 

The dependence of the kinetic friction force on tem- 
perature is shown in Fig. |31 For all "non-special" angles, 
when the static friction is relatively high, the kinetic fric- 
tion decreases with temperature (e.g., see the dependence 
for (p = 34.3° in Fig. HT), reflecting the standard ther- 
molubric effect [l|, |27H33| due to temperature-assisted 
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Figure 2: (Color online) Static (panel a, stick slip at « = 0.03) 
and kinetic (panel b, smooth sliding with v = 0.3) friction 
forces for several values of the misfit angle (j> for T — and 
different thicknesses of the lubricant film: Ni = 1 (diamonds), 
2 (green up triangles), 3 (blue down triangles), and 5 (red 
circles). Insets display the same dependences in logarithmic 
scale. All forces here and in the following figures are reported 
in the natural model units, i.e. in units of Vss/rss- 



barrier overcoming. However, for all special angles pro- 
ducing the SS regime, the behavior is different; thermal 
fluctuations perturb a rather delicate solitonic motion, 
leading to an initial increase of friction with T. Such a 
behavior occurs also for thicker films. For the misfit an- 
gle (j) ~ 7.3° we observe that at v = 0.1 the one-layer film 
reaches an exceptional sliding state characterized by very 
low friction due to SS; this type of superlubricity is not 
typical, and it is quickly destroyed with the increase of 
either temperature (Fig.|4]) or velocity (Fig. [3]). Note also 
that the superlubric SS regime is recovered significantly 
after it is abandoned as temperature is lowered (dotted 
line, open symbols), thus opening a nontrivial hysteretic 
loop in the thermal cycle. 

Consider now a thicker lubricant film with Ni = 2 (see 
Fig. [5]). For w > 1 we observe the LS regime, where the 
lubricant is 3D molten (except the "special" angles 4> = 
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Figure 3: (Color online) Friction force as a function of the 
misfit angle (j) for the one-layer lubricant film at T = 0. Dia- 
monds correspond to static friction, triangles to kinetic fric- 
tion at t; = 0.1, and circles to kinetic friction at t; = 0.3. 
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Figure 4: (Color online) Kinetic friction force as a function 
of temperature (in natural units of Vas/iSks), where /cb is 
the Boltzmann constant) computed at u = 0.1 for the mis- 
fit angles = 17.7° and = 40.5° (solitonic), (j> = 34.3° 
(high friction), and cf) = 7.3° (fragile solitonic) for the A'^; = 1 
system. 



25.7° and 40.5° at w = 1, where we have LoLS between 
the two attached lubricant layers). At lower velocities, 
V < 0.3, the behavior is as follows. For (/) = at stick, the 
two ordered lubricant layers are ordered and attached to 
the corresponding substrates, but are 3D melted at slips. 
For all other angles, the LoLS regime operates during 
slips: for (f) < 25° the attached layers are 2D molten, 
while for > 28° the attached layers remain ordered. 
Finally, for the angles (f> = 17.7°, 25.7° and 40.5° the 
friction forces exhibit deep minima produced by the SS 
solitonic mechanism. 

We come now to describe the results for the Ni — 3 
system as a prototypical lubricant of mesoscale thick- 
ness. The static friction, as well as the kinetic friction 
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Figure 5: (Color online) Friction forces versus the misfit angle 
(f> for the Ni = 2 system. Diamonds and solid curve describe 
the static friction, dotted curves show the kinetic friction force 
at different driving speeds, up to values comparable to the 
lubricant sound velocity. 



in the LoLS regime, can change by more than one order 
in magnitude when the misfit angle varies, as illustrated 
in Fig. m (j) — produces the highest friction like for 
thinner lubricants. The smooth sliding, as well as slips 
during stick-slip, correspond to either the LoLS or the 
LS regime. For all other angles 7^ 0, sliding always 
corresponds to the LoLS regime at low driving velocities 
V < 0.3. Contrary to the (j) = case, now sliding is typ- 
ically asymmetric and takes place at one interface only, 
between the middle layer and one of the attached layers, 
so that the middle lubricant layer sticks with either the 
top or the bottom substrate. The middle layer may re- 
main ordered during sliding or, for some values of 0, it is 
2D-melted; in the latter case, the friction force is higher. 
For the special misfit angles <j) = 17.7°, 25.7° and 40.5° 
identified by stars in Fig. [HI we again observe the "su- 
perlubricity" characterized by the very low friction. In 
these cases, the lubricant film remains solid and ordered 
during sliding, and moves as a whole with the top or bot- 
tom substrate, in a SS sliding. However, the lubricant is 
not rigid during motion, thus enhancing the "solitonic" 
mechanism. 

The results described above, remain qualitatively the 
same at nonzero temperatures T < Tm- For example. 
Fig. |6}d shows the dependence of friction force on for the 
"room" temperature T = 0.025. Both the static and ki- 
netic (for the LoLS regime) friction forces decrease when 
T increases. However, for the SS regime, the behavior is 
different - fiuctuations due to temperature perturb the 
solitonic motion, leading to an increase of friction. 

The thicker lubricant film, Ni — b, behaves similarly 
to A^; = 3, as illustrated in Fig. [T] The main difference 
is the lack of a maximum in /«((/') at = 0. Again, 
for the misfit angles = 17.7° and 25.7° we observe 
"superlubric" sliding. For smooth sliding with v < 0.3 as 
well for slips during stick-slip motion, we observe either 
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Figure 6; (Color online) Friction forces as functions of the 
misfit angle <jf> for the three-layer system. Diamonds and solid 
line represents static friction; stars and dotted lines mark the 
regime of solitonic "solid sliding" with very low friction ( "su- 
per lubricity"). Dashed lines show the kinetic friction force 
at different driving velocities: v — 0.1 (blue down triangles), 
0.3 (orange circles), 1 (cyan left triangles) and 3 (red right 
triangles). Panel (a) is for T — 0, and (b) is for the "room" 
temperature T = 0.025. 
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Figure 7: (Color online) Friction forces versus the misfit angle 
(/!> for the Ni — 5 system. Diamonds and solid curve describe 
the static friction, dotted curves show the kinetic friction force 
at driving velocities v = 0.1 (triangles) and 0.3 (circles). 



the LoLS regime, where the three central layers are 2D 
melted and sliding occurs between the middle layers (e.g., 
between layers 1-2 or 2-3), or the LS regime, where all 
three middle layers are 3D-molten. The LoLS regime 
marks the dips in f^, while LS generates larger /j., as 
occurs in the angular intervals 0° < (j) < 10° and 20° < 
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Figure 8: (Color online) Friction force versus the misfit angle 
(j> for the A^; = 1 system, driven at v = 0.01, for different 
values of the applied load 
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Figure 9: (Color online) The probability distribution of static 
thresholds for Ni = 1, 2, 3 and 5 (insets) and the averaged 
distribution Pc{fs)- These distributions are obtained in the 
assumption that all misfit angles (jy are equally likely. The 
calculated function Pc{fs) may be approximately fitted by 
the distribution Pdfs) = f'^fs exp{-fs/f) with / = 0.005. 



(j) < 25°. For larger velocities, f > 1, all five lubricant 
layers are melted and the LS regime operates in full. 

Figure [8] reports the friction force for 4 values of the 
applied load. These calculations demonstrate that the 
dependency of friction on the substrate rotation is very 
similar for different loads, with a general increase of fric- 
tion with load. We obtain very similar results for Ni ^ 2 
and 3. However, for thicker films, Ni = 5, the situa- 
tion changes dramatically at high load (// = 1): the film 
rearranges into a closely packed four-layer configuration 
which remains solid under sliding. As a consequence, the 
peak structure of friction as a function of cj) changes as 
well, because the lubricant structure acquires more atoms 
per layer and changes symmetry. 



Because the static friction varies over such a large in- 
terval, its specific value for a given misfit angle (j) has 
little importance. Indeed one could hardly control the 
misfit angle in a real system, except in especially favor- 
able situations [l^. Moreover, a system where the sliding 
surfaces have an ideal crystalline structure oriented with 
a controllable (j) is exceptional. Real surfaces usually have 
areas (domains) with different orientation. For polycrys- 
talline substrates, it is reasonable to assume that all an- 
gles are presented with equal likeliness. It may then be 
more interesting to examine the probability distribution 
of fs values, regardless of (f>. The insets of Fig. [3] report 
the histograms of forces as resulting from our simulations 
of different thicknesses. 

Besides, if we also assume that the lubricant film is not 
uniform but has different thickness at different places (it 
is certainly so if the surfaces have some roughness) , then 
we can average over different thicknesses; the resulting 
distribution Pdfs) is shown Fig. [9] It is precisely this 
distribution which represents the main output of our MD 
simulations, as it then allows us to predict tribological be- 
havior of the system with the help of a master-equation 
approach [33 | . The calculation summarized in Fig. IH] are 
carried out at // = 0.1. The effect of a load increase 
would be mainly to scale the Pdfs) distribution, so that 
it would peak at larger friction. The resulting distribu- 
tion displays several spikes which are likely due to the 
fixed size of the simulated contact. We will investigate 
the role of contact size on the Pdfs) distribution using 
a simplified model, in a separate publication [35j . 



IV. DISCUSSION AND CONCLUSIONS 

The main results of the present work can be summa- 
rized as follows: (i) The relative rotation of sliding con- 
tacts in a lubricated context promotes LoLS more fre- 
quently then the standard LS. (ii) For a few special an- 
gles LoLS leads to superlubric sliding completely anal- 
ogous of the unlubricated sliding of misaligned perfect 
crystalline substrates. This superlubric regime is however 
delicate and can be suppressed by small relative rotations 
of the substrates, by temperature, or by velocity-induced 
local heating, (iii) In a regime of boundary lubrication, 
friction forces do vary quite significantly with the rela- 
tive substrate relative orientation (f), even when the lu- 
bricant film becomes several atomic layers thick, (iv) To 
describe macroscopic friction in a context of multiasper- 
ity contact, where relative orientation is not really under 
control, the most important information to be extracted 
from MD simulations is a probability distribution Pdfs) 
rather than specific values of the static friction force fg. 

The present calculation is consistent with a rapidly 
(approximately exponentially) decaying distribution 
Pdfs), up to a cutoff force related to the average contact 
size. When temperature can promote rotations, beside 
the standard reduction of friction at low speed due to 
thermal crossing of barriers, thermal fluctuations could 




Figure 10: (Color online) Construction of the rotated lattice. 



affect the barrier distribution itself by suppressing small 
barriers in favor of higher ones j36| . 
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Appendix: The construction of the rotated substrate 

In this Appendix we report the construction of the 
substrate rotated by a given misfit angle (f>. The main 
problem here is to obtain numerous misfit angles (j) which 
satisfy PBC in x and y directions simultaneously, while 
maintaining a constant simulation size, and rectangu- 
lar PBC (for the square shape of the simulation cell 
the construction is much simpler, e.g., see Ref. [isj). 
The idea of the construction is demonstrated in Fig. (TU] 
The substrate is arranged according to a square lattice 
with lattice constant Ug, and the simulation cell area is 
X Ly — Mxtts X MyQs- The rotated bottom lattice 
is constructed as a set of parallelograms, so that the el- 
ementary cell of the new lattice has size x ay with a 
base angle Tr/2 + e. In the perfect case we would have 
and e = 0. However, to satisfy the PBC, 
the rotated lattice has to be distorted from the ideal 
square shape, and the idea is to reduce this distortion 
to a minimum. 

The rotated lattice is defined by two integers ni and 
712 (see Fig. fTU]) which determine the rotation angle (p. 





(a) (b) 



Figure 11; (Color online) Two typical examples of the rotated 
substrate lattice, with Mx = 12 and My — 11: (a) for (f) = 
28.75° (m = 10, 712 = 6, e = 0.58°, = 3.15, Uy = 2.86, 
AiVsub = -2, 5i = 0.0294), and (b) for <j) = 34.35° (ni = 
10, n2 = 7, £ = -2.3°, Ox = 3.01, ay = 2.99, AN.uh = 0, 
Si — 10~^). The atoms are shown by circles: those within the 
simulation cell are red/dark, while their periodic "images" are 
cyan/clear. 



For example, the choice ni = Mx and n2 = or ni = 
and n2 — My gives the original square lattice, while the 
sets Til = Mx and 712 = 1 or rii = 1 and 712 = My 
lead to the minimally allowed misfit angle for the given 
size of the simulation cell. Let us draw two lines (see 
Fig. [TU]), the first line starts at the point O2 = (0,ivy) 
and has the length niax, while the second line starts 
at the point O3 = {Lx,Ly) and has the length 7120^. 
These lines intersect at point O4 forming an angle Tr/2+e. 
The rotated substrate atoms are placed along these lines, 
and then periodic shifts by multiples of and ay in the 
directions defined by these two lines will generate the 
rotated oblique lattice. 

The oblique lattice constants Ux and ay are determined 
by two constrains. First, we must preserve the area of 
the elementary cell: 



axay cose 



Second, from the triangle O2O3O4 we have 

Lx — (niax)'^ + (772^1/)^ + 2771 a^ tane . 
From Eqs. (jA.ip and (|A.2p we obtain 



'1,2 



{A cos^ e±D) / {2nl cos^ e) 



(A.l) 



(A.2) 



(A.3) 



1 /2 

where D — (A^ cos^ e — Anfn'^a'^) cose and A = Lx ~ 
2nin2a^ tan e. The signs ± in Eq. (jA.3|) yield two possible 
solutions which we call as "left" and "right" ; typically we 
use the "right" variant when 771 > Mx/2 and the "left" 
one for 711 < Mx/2. Finally, the actual misfit angle (p is 
given hy (j) — If + e/2, where 



siniyj = 772 flj, cose/Lx 



(A.4) 



The construction described above guarantees the per- 
fect PBC in the x direction, but not in the y one. There- 
fore, a next step in construction is to characterize this 
distortion. Considering the triangle O1O3O5 in Fig. (TU] 
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(the line O1O5 is parallel to O2O4), the lengths of its 
short sides are — Rb sin((/j + e) and Ry = Rb cos ip, 
where Rb = Ly/ cose. The distortion of PBC is deter- 
mined by two parameters Sx = [Rx/cix — ^i^i{Rx / clx)]^ 
and Sy — [Ry/uy — iiit{Ry/ay)] . In the ideal case, it 
should be Sx ^ Sy ^ 0. 

A "quality" of the rotated lattice can be characterized 
by two parameters: the first parameter 

Si=Sl + Sl (A.5) 

describes the distortion of periodic boundary conditions, 
while the second parameter 

S2 = {ax/ay - 1)2 (A.6) 

indicates how close ax and ay are to the original square- 
lattice constant as . Perfect rotations are realized for Si ~ 



and ^2 = simultaneously at e = 0. Thus, for given 
integers ni and n2, we plot Si and S2 as functions of e, 
to choose an appropriate minimum of Si{e) close to the 
point £ = 0, and to check that S2 is not too large at that 
point. 

Note that in the non-ideal case, some atoms within the 
Lx X Ly simulation cell may be missing or, for other sets of 
parameters, some atoms may overlap with their "image" 
atoms generated by PBC. To overcome this problem, we 
shifted slightly the bottom boundary of the selected area. 
As a result, in the rotated lattice the number of substrate 
atoms may differ from the original one by a value AA'sub- 

Finally, because different sets of parameters may pro- 
vide approximately the same misfit angle, we can choose 
the best set, the one which minimizes Si^2 and e. Two 
typical examples of the construction described above are 
shown in Fig. [TT] 
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